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Abstract 

The sea slug Elysia chlorotica offers a unique opportunity to study the evolution of a novel function (photosynthesis) in a 
complex multicellular host. Elysia chlorotica harvests plastids (absent of nuclei) from its heterokont algal prey, Vaucheria 
Utorea. The "stolen" plastids are maintained for several months in cells of the digestive tract and are essential for animal 
development. The basis of long-term maintenance of photosynthesis in this sea slug was thought to be explained by 
extensive horizontal gene transfer (HGT) from the nucleus of the alga to the animal nucleus, followed by expression of 
algal genes in the gut to provide essential plastid-destined proteins. Early studies of target genes and proteins supported 
the HGT hypothesis, but more recent genome-wide data provide conflicting results. Here, we generated significant 
genome data from the E. chlorotica germ line (egg DNA) and from V. Utorea to test the HGT hypothesis. Our compre- 
hensive analyses fail to provide evidence for alga-derived HGT into the germ line of the sea slug. Polymerase chain 
reaction analyses of genomic DNA and cDNA from different individual £. chlorotica suggest, however, that algal nuclear 
genes (or gene fragments) are present in the adult slug. We suggest that these nucleic acids may derive from and/or reside 
in extrachromosomal DNAs that are made available to the animal through contact with the alga. These data resolve a 
long-standing issue and suggest that HGT is not the primary reason underlying long-term maintenance of photosynthesis 
in E. chlorotica. Therefore, sea slug photosynthesis is sustained in as yet unexplained ways that do not appear to endanger 
the animal germ line through the introduction of dozens of foreign genes. 
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Introduction 

Many species of herbivorous sacoglossan molluscs (sea slugs) 
exhibit Septoplasty, which is the uptake and retention of 
plastids as a result of feeding on algal prey (see reviews by 
Trench 1975; Jensen 1997; Yamamoto et al. 2009; Rumpho 
et al. 2011; Wagele et al. 2011; Pierce and Curtis 2012). 
Individual kleptoplastic associations vary widely in terms of 
specificity of the animal toward its algal prey and in retention 
time and functionality of the captured plastids (see reviews 
by Evertsen et al. 2007; Evertsen and Johnsen 2009; Handeler 
et al. 2009; Wagele et al. 2011; Rumpho et al. 2011; Maeda 
et al. 2012). Elysia chlorotica is at one extreme because its 
known algal prey is limited to two marine species of the 
coenocytic heterokont genus Vaucheria (i.e., Vaucheria Utorea 
and V. compacta; West 1979; unpublished observations), and 
because it retains functional, nondividing plastids in the cells 
lining the digestive tract for several months (Pierce et al. 1996; 
Green et al. 2000). The long- term functioning of V. Utorea 
plastids in £ chlorotica is perplexing because no evidence 



exists for the presence of algal nuclei in the sea slug, and 
the majority of plastid proteins are encoded in the nuclear 
genome, in addition to all the nuclear regulatory factors that 
are targeted to the plastid (Rumpho et al. 2008; Woodson and 
Chory 2008; Puthiyaveetil et al. 2010). 

Several hypotheses have been proposed to explain how the 
plastids continue to function within the animal cells, includ- 
ing horizontal gene transfer (HGT) of essential photosynthesis 
genes from the nuclear genome of the algal prey to the sea 
slug nuclear genome. The animal digestive cells presumably 
come into contact with algal nuclear DNA during suctorial 
feeding on Vaucheria, although mechanisms to explain trans- 
fer of HGT-derived genes to the germ line of the animal 
have not yet been proposed. The HGT hypothesis in E. chlor- 
otica has been driven by biochemical data (Pierce et al. 1996, 
2010; Rumpho et al. 2008, 2009; Soule and Rumpho 
2012), polymerase chain reaction (PCR), and real time 
quantitative reverse transcription-PCR (qRT-PCR) data 
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(Pierce et al. 2007; Rumpho et al. 2008, 2009; Schwartz et al. 
2010; Soule and Rumpho 2012; present study), and most 
recently by the analysis of £ chlorotica transcriptome data 
(Pierce et al. 2012). 

Sequence analysis of transcript data from multiple adult 
photosynthetic £ chlorotica has provided conflicting evi- 
dence for Vaucheria- or other algal-derived transcripts in 
this animal (Pelletreau et al. 2011; Pierce et al. 2012). In the 
initial transcriptome studies carried out by our laboratory, 
no evidence was found in the £ chlorotica raw or assem- 
bled transcripts for horizontally transferred algal genes 
(Pelletreau et al. 2011). This was based on Basic Local 
Alignment Search Tool (BLAST) analysis against the nonre- 
dundant National Center for Biotechnology Information 
(NCBI) database and also by using Phylosort, a phyloge- 
nomic pipeline designed to identify candidates of gene 
transfer (Moustafa et al. 2009). Although a handful of 
putatively transferred genes (mostly of bacterial or viral 
origin) emerged from these analyses, they were neither 
algal derived nor involved in photosynthetic processes. 
This study complemented the results obtained by Wagele 
et al. (2011) whose partial transcriptome data from two, 
albeit very different kleptoplastic sea slugs (£ timida and 
Plakobranchus ocellatus), also revealed no evidence for 
HGT. In contrast, Pierce et al. (2012) reported the presence 
of 111 algal transcripts, 52 of which are the products of 
nuclear genes, in the adult nuclear genome of £ chlorotica. 
Although each of these transcripts was only found as a 
singleton read in the animal transcriptome library, and 
none were greater than 90 bp in length, these authors pre- 
sumed the presence of these transcripts indicated HGT 
from the algal nucleus and integration into sea slug geno- 
mic DNA (Pierce et al. 2012). More recent data clearly 
indicate the presence of algal transcripts in the form of 
singletons in the animal transcriptome; however, no data 
to date provide strong evidence for integration of intact 
algal nuclear genes into the sea slug genome or, more 
specifically, into the germ line (also see discussion with 
supplementary fig. S1, Supplementary Material online). 

Given these conflicting data, we took the direct 
approach of generating significant genome data from nat- 
urally aposymbiotic £ chlorotica egg DNA that had never 
come in contact with algal prey. We analyzed individual 
genome reads using BLAST methods to search for evidence 
of HGT. To strengthen our analysis, we also generated 
transcriptome data from the algal prey V. litorea and 
used these sequences to query the egg DNA to look for 
instances of DNA transfer. Finally, we used all the putative 
HGT transcript data from Pierce et al. (2012) as a query 
against egg DNA. Taken together, the results of our anal- 
yses indicate the absence of algal-derived genes in the germ 
line of £ chlorotica. Therefore, all these results point to yet- 
to-be-explained mechanisms for the long-term mainte- 
nance of photosynthesis in the gut of this animal. Here, 
we present initial results from the PCR analysis of individual 
animals that suggest the existence of a novel mechanism of 
transient gene expression in adult sea slugs that may help 
explain the conflicting data. 



Results and Discussion 

Looking for Hits to V. litorea Transcripts in the 
E. chlorotica Genome 

Sequencing of the V. litorea transcriptome on an lllumina 
MiSeq instrument generated 21.3 million paired-end reads 
(average size=151bp) totaling 3.2 Gbp. These data were 
assembled into 65,468 contigs (N50 = 824bp, median cover- 
age = 2.3 x, total assembly size = 25.5 Mbp). The contigs were 
analyzed using BLASTx (e-value cutoff < 10 -5 ) against our 
comprehensive local database (e.g., Chan et al. 2011; Price 
et al. 2012). This database consists of 16.9 million protein 
sequences derived from RefSeq v.51 with the addition of 
sequenced eukaryote (e.g., fungal, metazoan, viridiplant, and 
stramenopile) genomes from the Joint Genome Institute 
(http://www.jgi.doe.gov/, last accessed May 16, 2013) and 
six-frame translated eukaryote expressed sequence tag (EST) 
sequences obtained from the NCBI dbEST (http://www.ncbi. 
nlm.nih.gov/dbEST/, last accessed May 16, 2013). The data 
showed that 17,311 EST contigs encode proteins with 
hits to the database with many showing the 
expected "chromalveolate" (e.g., diatoms and oomycetes 
[stramenopiles]) affiliation (fig. 1; supplementary table S1, 
Supplementary Material online). We also generated 153.2 mil- 
lion paired-end lllumina reads of average size 146 bp from 
£ chlorotica total egg genomic DNA, totaling 22.5 Gbp of 
data. BLASTn analysis using the assembled V. litorea cDNA 
sequences as a query against the sea slug egg genome database 
at an e-value cutoff of <10 -5 yielded 1,298 top hits to 241 
nonredundant reads from the sea slug (supplementary table 

52, Supplementary Material online). Three sea slug genomic 
sequences that were inserted into the algal query data set as 
positive controls were returned in the search, indicating that 
our BLAST parameters were appropriate for a large database 
of short-read data (results not shown). 

Although a large number of hits were returned between 
the algal EST contigs and animal singletons, the average 
length of these hits was 36 bp (four were > 100 bp), 
and only 76 had an e-value less than 10 -20 . Nonetheless, to 
investigate the function of these algal hits, each sequence 
was annotated using the NCBI nr database. This analysis 
showed that 730 had no significant hit to the NCBI database, 
and the remaining 568 contigs were comprised predomi- 
nantly of housekeeping and/or conserved genes (e.g., actins, 
tubulins, and some hypothetical proteins; see supplementary 
table S2, Supplementary Material online; fig. 2). Although 
this comparison generated numerous top hits to strameno- 
piles, this result reflects the phylogenetic affiliation of the 
V. litorea query sequence and does not indicate that the 
homolog in £ chlorotica is of algal origin. Indeed, when 
the corresponding £ chlorotica singletons (241 were unique 
of the 1,298 matches) were blasted against the NCBI nr data- 
base, 96 were annotated and many of the top hits were 
to metazoan or bacterial sequences (supplementary table 

53, Supplementary Material online; fig. 2). A handful of 
genes appeared to be of plant origin; however, on closer 
inspection, each of these genes aligned to the same 
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Number of Hits 




Stramenopiles 
Metazoa 
Viridiplantae 
Fungi 
Alveolata 
Excavata 
Amoebozoa 
Cryptophyta 
Proteobacteria 
Haptophyceae 
Rhizaria 
Rhodophyta 
Choanoflagellida 
Cyanobacteria 
Firmicutes 
Bacteriodetes 



Fig. 1. BLASTx results using assembled Vaucheria litorea cDNA contigs 
to query a comprehensive local genome database. This database consists 
of 16.9 million protein sequences derived from RefSeq v.51 plus 
sequenced eukaryote (e.g., fungal, metazoan, viridiplant, and 
stramenopile) genomes from the Joint Genome Institute (http://www. 
jgi.doe.gov/, last accessed May 16, 2013) and six-frame translated eu- 
karyote EST sequences obtained from the NCBI dbEST (http://www. 
ncbi.nlm.nih.gov/dbEST/, last accessed May 16, 2013). 



singleton (see supplementary table S3, Supplementary 
Material online; HWUSI-EAS1 588:2:1 0:1 3875:1 51 05#0/3) de- 
rived from £ chlorotica, which, alone, had no significant sim- 
ilarity to any gene in NCBI. Therefore, in spite of the finding of 
regions of limited similarity between the algal transcriptome 
and the animal genome (as might be expected for such com- 
plex data), we conclude that genomic DNA from £ chlorotica 
eggs is devoid of any full-length (or identifiable, partial) genie 
regions from V. litorea. Most importantly, our analysis does 
not return any genes or DNA regions shared by V. litorea and 
the sea slug that are involved in photosynthesis. 

Examining Possible Instances of HGT in the 
£ chlorotica Genome 

The results described earlier, combined with our previous 
work (Pelletreau et al. 2011), stand in stark contrast to the 
results of Pierce et al. (2012). These latter authors compared a 
transcriptome from pooled adult, starved £ chlorotica against 
a draft assembly of the genome of V. litorea and with EST 
contigs (independently derived from our data) from this alga, 
and identified 111 HGT candidates (derived from 431 ESTs). 
Of these 111 animal-derived individual transcripts (out of a 
total of 98,238,204 EST reads [i.e., 0.00011% of the total read 
count]), 52 were nuclear-encoded plastid proteins with 27 
involved directly in photosynthesis. Given these results, we 
searched the 22.5 Gbp of raw sea slug egg genome reads 
that we generated for hits to any of the existing photosynthe- 
sis genes that have been identified in previous transcriptome 
analyses of adult £ chlorotica ("HGT Library" fig. 2C). These 
included the 431 transcripts identified by Pierce et al. (2012) 
and the psbO and prk transcripts identified using targeted 
PCR (Rumpho et al. 2008, 2009; Schwartz et al. 2010). 
We found BLASTn hits (e-value cut off of <10~ 5 ) for 177 of 
the 431 transcript sequences from Pierce et al. (2012) in our 



£ chlorotica egg genome library (supplementary table S4, 
Supplementary Material online); however, none of the hits, 
aside from an HSP70 molecular chaperone, had a blast 
score more than 59 or an alignment length of greater 
than 46 nucleotides (see example in supplementary fig. S2, 
Supplementary Material online). All three positive con- 
trol sequences were returned in the BLAST output. Thus, 
we were unable to identify any of the putative algal-derived 
examples of HGT in individual reads from the egg genome. 

Testing the Robustness of Our Analysis 
One commonly cited reason for the lack of evidence of HGT 
in transcriptomic and genomic analysis is insufficient cover- 
age resulting in the failure to detect rare transcripts or low 
copy number genes (e.g., Pierce et al. 2012). To test the com- 
pleteness of our egg genome sequence data, we mapped the 
reads onto the £ chlorotica mitochondrial genome (Rumpho 
et al. 2008) using the CLC Genomics Workbench (similar- 
ity =0.90, length fraction = 0.90). This analysis showed that 
the egg data spanned 14,132 bp (100%) of the organelle 
genome with an average coverage of 1,976 x. We repeated 
the analysis with the nuclear-encoded histone H3 and actin 
sequences available in Gen Bank (accession nos. JN819183 and 
AF448493) and found our data span 100% of each coding 
region at average coverage of 10,319x and 126x, respectively. 
Finally, we mapped the £ chlorotica egg genome reads against 
a subset of the assembled EST contigs from Pelletreau et al. 
(2011) that had a length >300 bp (i.e., a total of 16,012 con- 
tigs). This analysis showed that 16.01 million reads were 
aligned (at a stringency of 92.5% similarity over 55% coverage; 
the latter to allow for introns in the egg genome) to 15,771 
(98.5%) of the EST contigs. The median coverage of the map- 
ping was 19.23x, and the average coverage was 198x. These 
results suggest that the majority of the animal gene inventory 
is present in our egg genome data and, therefore, is present in 
the germ line. The coding regions we have identified occurred 
in gene-rich regions that have more than 100x average 
genome read coverage. Nevertheless, we stress that our results 
do not prove the absence of HGT candidates in the nuclear 
genome but do, however, argue strongly against their 
presence. To unambiguously prove this point will require a 
completed sea slug genome sequence. 

Previous studies have demonstrated the presence of plas- 
tid-encoded gene transcripts in adult sea slugs, even though 
polyadenylated transcripts were targeted to exclude much of 
the plastid-derived RNA during library construction (Mujer 
et al. 1996; Pelletreau et al. 2011; Pierce et al. 2012). 
The presence of plastid transcripts is nevertheless expected 
because sea slugs maintain photosynthesis for protracted 
periods of time (at least 6 months). Given this observation, 
it was also crucial to show that our egg DNA assembly did not 
contain plastid genes. To address this issue, we used BLASTn 
and the completed V. litorea plastid genome (Rumpho et al. 
2008) as a query against the egg DNA data. This analysis 
did not turn up a single significant hit (e-value cut off 
<10 -5 ), indicating the absence of algal plastid DNA (i.e., 
lack of algal contamination) in our sea slug egg genomic 
DNA preparation. 
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Fig. 2. Comparative gernomics approach used to test for the presence of alga-derived genes (i.e., HGT candidates) in the germ line of Elysia chlorotica. 
(A) Algal transcriptome generated in this study. (B) Sea slug egg genomic data generated in this study. (C) Cummulative library of HGT candidates from 
past studies. Results of the various comparisons are illustrated and listed in detail in the indicated supplementary tables, Supplementary Material online. 
No algal or photosynthetic genes were detected in any of the queries performed with the £. chlorotica genomic reads. 



Explaining the Inconsistency in Evidence for 
HGT in Sea Slugs 

As discussed earlier, previous work has suggested in E. chlor- 
otica the presence of genomic DNA or cDNA encoding 
Vaucheria-derived genes involved in photosynthesis, not 
only in adult animals but also in veliger and egg DNA (for 
discussion, see Pierce et al. 2012). Given that our analysis of 
genomic and EST data do not provide any support for the 
presence of algal genes of any kind in the germ line of the sea 
slug it becomes important to explain the findings of algal 
genes in previous studies. The absence of algal genes found in 
our analysis of the egg genome and in a comprehensive anal- 
ysis of the V. litorea cDNAs leaves at least two possible expla- 
nations for the inconsistency in current data. The first is that 
contamination of algal DNA or cDNAs in the PCR reactions 
may have provided artifactual evidence for Vaucheria se- 
quences in the sea slug This explanation seems on the surface 



to be unlikely because the PCR experiments have been done 
by different investigators over the course of several years and 
have relied on different DNA/cDNA sources and PCR re- 
agents. The second explanation is that the trace amounts 
of algal DNA recovered using PCR may indicate true HGT 
events that have simply gone undetected in our draft assem- 
bly. If, however, our egg genome data are an accurate repre- 
sentation of this compartment, then it is also possible that the 
algal DNAs are not integrated into the egg genome but rather 
may reside in extrachromosomal or other forms of DNA that 
were not sequenced using our approach (see Discussion 
later). Therefore, regardless of whether evidence is ultimately 
found for algal DNAs in the sea slug genome, the overwhelm- 
ing evidence until now suggests that, consistent with a grow- 
ing list of other studies, no evidence exists for massive HGT 
from partner symbionts into chromosomal DNA of the host 
(e.g., Wisecaver and Hackett 2010; Pelletreau et al. 2011; 
Wagele et al. 2011; Pillet and Pawlowski 2013). An interesting 
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exception is the recent finding of gene acquisition "from the 
environment" by bdelloid rotifers that provides this lineage 
with novel physiological adaptations (Boschetti et al. 2012). 

A Final Note on Study Design and PCR Data 
We suggest that the design of existing studies may help to 
explain the conflicting data regarding alga-derived HGT in 
the sea slug. In particular, previous studies have relied on 
the use of multiple pooled animals rather than individuals 
when characterizing HGT. Here, using PCR primers that target 
different HGT candidates (table 1), we found an inconsistent 
pattern with respect to the presence or absence of candidate 
genes between different individuals (fig. 3) and also over 
time (fig. 4). Using these preliminary data, we identified 
three major trends: 1) there is high variability in the pres- 
ence/absence of HGT-implicated genes among individual 
animals, regardless of treatment (i.e., time spent feeding or 
breeding history), 2) the presence of a gene is not always 
consistent with the presence of transcripts from the same 
individual, supporting very low copy number of any particular 
transcript, and 3) traditional PCR continues to provide evi- 
dence that these genes are present in photosynthesizing 
adults, suggesting their role in individual animal growth and 
development. By using individuals and looking at exposure to 
algal DNA as variables, we identified a dynamic process of 
gene presence and expression existing in individual adult an- 
imals. For example, if one were to pool the three biological 
replicates used in figure 4, then all 11 of these genes would be 
"present" in the DNA leaving the false impression that HGT is 
present in all members of the sea slug population. In earlier 
work with this system, PCR-based approaches provided pos- 
itive results for only 15-20% of the samples (Pierce et al. 
2007), which may be explained by this large amount of in- 
terindividual variability. 

A possible explanation for the unstable pattern of "HGT" 
detected here is, as suggested earlier, the presence of extra- 
chromosomal DNA in £ chlorotica that may encode foreign 
genes and/or gene fragments. Extrachromosomal DNA has 



been reported in a number of organisms ranging from pro- 
tozoa to vertebrates (Cohen et al. 2007; Cohen and Segal 
2009), and in several cases, its formation and/or disappear- 
ance is linked to specific developmental changes (Noto et al. 
2003), and more recently, to distinct gene regions (Shibata 
et al. 2012). It is, therefore, possible that the acquisition of 
DNA from algal food in the sea slug differs among individ- 
uals each generation and is correlated with feeding. The gut 
of adult feeding £ chlorotica may function as an in vivo 
transfection system, whereby gene uptake is variable and 
random, often transient and extrachromasomal. In 
Caenorhabditis elegans, the model holocentric organism, 
extrachromasomal DNA can be heritable once integrated 
(albeit into the germ line; Stinchcomb et al. 1985), often 
in a mosaic pattern and propagated to daughter cells 
(Frokjaer-Jensen et al. 2008). Recent data suggest that nolo 
centrism may be more widespread than previously thought 
(Melters et al. 2012). 

Our work with the egg genome data suggests that foreign 
genes, either integrated or as extrachromosomal DNA, were 
not present in these data. In the latter case, the total DNA 
preparation we used for the egg genome would be expected 
to contain extrachromosomal material, including mini-cir- 
cles (i.e., gradient fractionation was not included during 
DNA isolation). This idea is supported by the finding of 
many copies of mitochondrial and bacterial DNA in the 
egg genome data. However, the copy number of extrachro- 
mosomal DNA could be below detection limits compared 
with the genomic (or mitochondrial and bacterial) DNA. 
The DNA used to generate the £ chlorotica genome data 
was isolated from multiple egg ribbons produced by multi- 
ple pairs of animals, resulting in thousands of individuals 
contributing their nucleic acids, and yet we did not observe 
any evidence of HGT. Although we cannot rule out that all 
these egg ribbons were unique in not possessing any HGT 
candidates, we can conclude that HGT is not universal 
throughout the population of £ chlorotica collected from 
Martha's Vineyard, MA. 



Table 1. Eleven Candidate Nuclear-Encoded Algal Genes That Have Been Implicated in HGT in Elysia chlorotica and Two Positive Control Genes 
(one algal and one animal) Used to Assess Interindividual Variability and Change over Time. 



Gene Target 


Primer Name 


Primer Sequence Source 


Function 


Spermadine synthase 


Spds 


Pierce et al. (2009) 


Algal nuclear control 


Beta actin 


P-actin 


Present study 


Animal nuclear control 


Light harvesting protein 1 


Lhd 


Pierce et al. (2010) 


Light harvesting 


Light harvesting protein 2 


Lhc2 


Pierce et al. (2010) 


Light harvesting 


Light harvesting protein 3 


Lhc3 


Pierce et al. (2010) 


Light harvesting 


Light harvesting protein 4 


Lhc4 


Pierce et al. (2010) 


Light harvesting 


Fucoxanthin binding protein 


fcp 


Pierce et al. (2007) 


Light harvesting 


Uroporphyrinogen decarboxylase 


uroD 


Pierce et al. (2010) 


Chlorophyll synthesis 


Chlorophyll synthetase D 


chlD 


Pierce et al. (2010) 


Chlorophyll synthesis 


Chlorophyll synthetase G 


chIG 


Pierce et al. (2010) 


Chlorophyll synthesis 


Chlorophyll synthetase H 


chlH 


Pierce et al. (2010) 


Chlorophyll synthesis 


Manganese stabilizing protein 


psbO 


Rumpho et al. (2008) 


Oxygen evolution 


Phosphoribulokinase 


prk.S2 


Soule (2009) 


Carbon fixation cycle 



Note. — All primers were targeted to exon regions to successfully amplify both genomic DNA and cDNA of the corresponding gene. 
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A WILD B LAB-REARED 




Fig. 3. PCR analysis showing the presence or absence of genes (A and B) or transcripts (C and D) of commonly identified HGT candidates in individual 
Elysia chlorotica collected from the wild (A and C) or 6-month-old laboratory-reared animals (B and D). Identical colored boxes represent the same 
individual being tested in panels A and C, and in panels B and D. Targeted genes are listed in table 1. 



In summary, despite the apparent complexity of the sea 
slug system, our data answer the fundamental question in this 
field: Is the ability of sea slugs to maintain algal plastids ex- 
plained by massive HGT to its germ line nucleus? Our work 
convincingly demonstrates that this is not the case. Elysia 
chlorotica does not appear to "endanger" its germ line by 
introducing dozens of foreign algal genes to support its pho- 
tosynthetic lifestyle. This leaves the far more interesting pos- 
sibility that the animal is able to harvest plastids and carry out 
photosynthesis using nonchromosomally integrated algal-de- 
rived genes, modified animal genes, or some combination of 
the two. Exploring these mechanisms using a combination of 
developmental ly targeted mRNA-seq approaches, proteo- 
mics, and deeper genome sequencing offer potential avenues 
to elucidate the (nearly abominable) mystery of long-term sea 
slug photosynthesis. 

Materials and Methods 

Collection of Sea Slug Tissue and Preparation of 
Egg Genomic DNA 

Specimens of £ chlorotica were collected from Martha's 
Vineyard, MA, in November 2009. Animal pairs were isolated 
and maintained at room temperature in the absence of algal 
prey in 0.2 Jim filtered, autoclaved 32 practical salinity units 
artificial sea water (ASW; 925 mosmol kg -1 Instant Ocean, 



Aquarium Systems, OH) and observed for mating. The sub- 
sequent egg ribbons were removed immediately upon laying, 
dipped in a 5% (v/v) sodium hypochlorite solution and rinsed 
in fresh water. The water was aspirated, and the eggs were 
flash frozen in liquid nitrogen. DNA was extracted using the 
Qiagen Plant DNeasy Kit (Qiagen Inc., Valencia, CA) and con- 
centrated by precipitation with ammonium acetate. Illumina 
paired-end sequencing libraries were prepared using the 
Nextera sample prep kit (Illumina Inc., San Diego, CA) before 
sequencing on an Illumina Genome Analyzer llx instrument 
using TruSeq SBS v5-GA reagents (Illumina Inc., San Diego, 
CA) following the manufacturer's protocols. Two flowcell 
lanes were sequenced in a paired-end (150 bp x 150 bp) 
manner. 

Algal Culturing and Preparation of RNA 
RNA was extracted from V. litorea maintained at room tem- 
perature at + 5 h into the light cycle (~30 u£/m 2 /s) using the 
Qiagen Plant RNeasy Kit. Illumina paired-end sequencing li- 
braries were prepared using the TruSeq RNA Sample 
Preparation Kit v2 (Illumina, Inc., San Diego, CA) before se- 
quencing on an Illumina MiSeq instrument using the MiSeq 
Reagent Kit (Illumina, Inc., San Diego, CA) in a paired-end 
(150 bp x 150 bp) approach. 
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Controls OHPF 1HPF 24HPF 3MS 



gDNA 



Controls OHPF 1HPF 24HPF 3MS 

cDNA 



Fig. 4. Percentage of animals (fl = 3) showing the presence of genes (left) or transcripts (right) encoding 11 putative HGT candidates (divided into 
functional categories) and two control genes (actin and spermadine synthase). Tested animals were starved for 3 months (3MS) and then allowed to 
feed on Vaucheria litorea for 1 h (OHPF) and subsequently starved for 1 h (1HPF) or 24 h (24HPF). Controls indicate samples that had combined algal 
and animal nucleic acids as templates. Targeted genes are listed in table 1. 
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Bioinformatic Analysis of Genomic and 
Transcriptomic Data 

The E. chlorotica egg DNA and V. litorea RNAseq reads were 
quality and adapter trimmed using the CLC Genomics 
Workbench (CLC Bio, Aarhus, Denmark). Elysia chlorotica 
read data were then exported in fasta format and converted 
to a BLAST database (Altschul et al. 1990). The V. litorea data 
were assembled into cDNA contigs using the Genomics 
Workbench and then "seeded" with three single read se- 
quences from the sea slug genome sequencing to act as pos- 
itive controls: one 146-bp sequence of 100% similarity, one 
truncated 75-bp sequence of 100% similarity, and one 75-bp 
sequence of 90% similarity. The algal contigs were then di- 
vided into 640 subsets and sent to a 640-core Beowulf style 
parallel computing cluster operated by the Rutgers University 
School of Environmental and Biological Sciences where each 
was used as a BLASTn (Altschul et al. 1990) query to the sea 
slug read database (e-value cutoff < 10 _1 ). An additional 
analysis was performed using the 431 individual plastid and 
nuclear-encoded transcript sequences of Pierce et al. (2012), 
some of which were hypothesized to be genes horizontally 
transferred from V. litorea to the E. chlorotica nuclear genome. 
These sequences were used as a BLASTn query (e-value 
cutoff < 10 _1 ) to the E. chlorotica egg genome database cre- 
ated above after insertion of three positive control sequences 
taken from the Pierce et al. (2012) data set: one 90-bp 
sequence of 100% similarity, one truncated 50-bp se- 
quence of 100% similarity, and one 90-bp sequence of 80% 
similarity. 

PCR Analysis of Individual E. chlorotica Specimens 
Interindividual Variability Experiments 
Wild animals were collected from Martha's Vineyard and 
maintained in the laboratory for more than 3 months without 
feeding. Laboratory-reared animals were raised as described 
by Pelletreau et al. (2012) and maintained at 10 °C on a 12:12 
L:D cycle being fed V. litorea weekly. 

Feeding Experiment 

Wild-caught animals maintained unfed in the laboratory 
for 3 months were placed in individual Petri dishes filled 
with ASW at room temperature on a 12:12 L:D cycle and 
allowed to acclimate for 48 h. The experiment was then 
initiated 45min after illumination by placing each animal 
on 0.15 g wet weight of V. litorea. Individual animals ob- 
served to feed continuously for 1 h were removed from 
V. litorea, rinsed in fresh water, blotted dry, and frozen 
in liquid nitrogen immediately after feeding (0HPF) or 
moved to new ASW and washed and frozen 1 h or 24 h 
after feeding ceased (1HPF and 24HPF, respectively). One 
set of animals was frozen before exposure to food and 
represented the state of animals after 3 months of starva- 
tion. Samples were stored at — 80 °C until further 
processing. 

DNA Extraction 

Total DNA was extracted from frozen V. litorea samples using 
DNAzol reagent (Invitrogen Co., Carlsbad, CA). Elysia 



chlorotica DNA was recovered during TRIzol RNA isolation 
(below) by following the published TRIzol DNA precipitation 
protocol to enable isolation of DNA and RNA from the same 
individual. 

RNA Isolation 

Total RNA was extracted from frozen algal samples using 
the Qiagen RNeasy Plant Mini Kit (Qiagen Sciences, 
Germantown, MD) following the manufacturer's protocol 
with all centrifugation steps performed at 4°C and 
8,000 x g. Total RNA was isolated from frozen sea slug 
samples using a modified TRIzol (Sigma, St. Louis, MO) 
procedure (Soule 2009). 

cDNA Generation 

High-quality RNA (0.5 jig) was treated with DNase I 
(Amplification Grade; Invitrogen), analyzed in an agarose 
gel to confirm the absence of genomic DNA, and first- 
strand cDNA was generated using Invitrogen's oligo-dT 
primed Superscript III First Stand Synthesis System. 

PCR Amplification 

PCR reactions were set up for 11 different algal nuclear- 
encoded plastid targeted genes and two control genes 
(table 1). All PCR reactions (except for spermadine 
synthase) were performed using 10 ng of template and 
1.25 U Promega GoTaq enzyme, 1x GoTaq Flexi buffer, 
0.4 jiM primers, 1.5 mM MgCI 2 , and 0.2 mM each dNTP, 
per reaction. Spermidine synthase required a final MgCI 2 
concentration of 3 mM. PCR conditions were optimized as 
follows: initial denaturing (95 °C, 2min), 40 cycles of 
touchdown (95 °C, 30 s, 60-50 °C for 1 min with a de- 
crease of 0.5 °C per cycle, and 72 °C for 45 s), and a 
final elongation step (72 °C, 10 min). PCR conditions for 
E. chlorotica /3-actin and Ihcal were only run on a 30- 
cycle touchdown program due to highly robust amplifica- 
tion. For every gene tested, biological replicates (n = 3) 
were used in each PCR along with no reverse transcrip- 
tase (no RT) controls (cDNA) and negative controls that 
lacked template DNA. 

Validation of PCR Results 

Products were analyzed with gel electrophoresis and 
sequenced by the University of Maine DNA Sequencing 
Facility to validate observed products. Sequences were com- 
pared against the NCBI nucleotide and protein databases 
using a combination of BLASTn and BLASTx, respectively 
(Altschul et al. 1990). 

Enhanced PCR to Verify Negative Trends 
For each gene targeted, cDNA reactions and their corre- 
sponding no RT controls for one replicate were reamplified 
for an additional 20 cycles to determine whether seemingly 
absent transcripts were present but at levels undetectable 
by standard PCR and gel electrophoresis/EtBr staining. The 
previously described touchdown PCR was run on these tem- 
plates for each time point as described earlier. Subsequently, 
1 jil of this PCR reaction was used as a template for an 
additional 20 cycles of reamplifkation at 50 °C annealing 
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temperature. PCR products were not observed in any of the 
reamplifkations. 

Supplementary Material 

Supplementary tables S1-S4 and figures S1 and S2 are avail- 
able at Molecular Biology and Evolution online (http://www. 
mbe.oxfordjournals.org/). 
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